Sufficient conditions for thermal rectification in graded materials 
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We address a fundamental problem for the advance of phononics: the search of a feasible thermal 
diode. We establish sufficient conditions for the existence of thermal rectification in general graded 
materials. By starting from simple assumptions satisfied by the usual anharmonic models that 
describe heat conduction in solids, we derive an expression for the rectification. The analytical 
formula shows how to increase the rectification, and the conditions to avoid its decay with the 
system size, a problem present in the recurrent model of diodes given by the sequential coupling 
of two or three different parts. Moreover, for these graded systems, we show that the regimes 
of non-decaying rectification and of normal conductivity do not overlap. Our results indicate the 
graded systems as optimal materials for a thermal diode, the basic component of several devices of 
phononics. 
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INTRODUCTION 

The study of the macroscopic laws of thermodynamic 
transport from the underlying microscopic models is still 
a challenge in statistical physics. In particular, the in- 
vestigation and control of the energy transport, which 
mainly involves conduction of heat or electricity, is a 
fundamental problem of huge theoretical and practical 
interest. The invention of transistor and other devices 
used to control the electric charge flow has led to the well 
known development of modern electronics. Its much less 
developed counterpart - the study and control of heat 
current - has, recently, presented interesting progress, 
promising to establish, in addition to electronics, a new 
physical branch in energy and information processing - 
the phononics [l|, Q : researchers have proposed nanode- 
vices such as thermal diodes or rectifiers 0-0] (already 
built in practice 6]), thermal transistors 0, thermal logic 
gates [8| and memories ^9]. The most fundamental com- 
ponent of these instruments is the thermal diode, a de- 
vice in which heat flows preferably in one direction. In a 
short analysis, we may say that this promising advance 
of phononics is directly dependent on the development 
of its basic component: a thermal diode with suitable 
properties. 

There are analytical attempts to investigate the phe- 
nomenon of thermal rectiflcation such as the works on 
spin-boson junctions billiard systems 11 1, etc., but 
most of the results are by means of computer simula- 
tions, see e.g. the work of B. Li and collaborators 12 1. 
The most common and recurrent design of diodes is given 
by the sequential coupling of two or more chains with dif- 
ferent anharmonic potentials Although frequently 
studied, this procedure is criticized due to the dif- 
flculty to construct such diode in practice, and due to 
the significative decay of the rectification with the sys- 
tem size. Recently, a different procedure was considered 
by Chang et al. |6|, who built, in an experimental work. 



the first microscopic solid-state thermal rectifier by using 
a graded material: nanotubes externally and inhomoge- 
neously mass-loaded with heavy molecules. It is worth 
to recall that graded materials, i.e., inhomogeneous sys- 
tems whose composition and/or structure change gradu- 
ally in space, are abundant in nature, can also be man- 
ufactured, and have attracted great interest in many ar- 
eas [l3|: there are many works devoted to the study of 
the electric, optical, mechanical and other properties of 
graded materials, but there are few studies in relation to 
their heat conduction investigation. 

In the present work, we address this fundamental prob- 
lem of phononics: the built of an appropriate thermal 
diode, namely, a simple system that may be constructed 
in practice, and with a rectification that does not decay 
with the system size. We start from simple conditions 
for the local thermal conductivity, conditions that are 
quite general and that are satisfied by anharmonic crys- 
tal models used to describe heat conduction in solids, and 
then we show that they are sufficient to lead to rectifi- 
cation in graded models. Moreover, we derive an expres- 
sion for such rectification that allows us to see how to 
make it larger, and how to avoid its decay with the sys- 
tem size. In short, we show that properly manipulated 
graded materials have suitable properties of rectiflcation, 
and so, they shall play a central role in the building of 
thermal nano-devices. The simplicity of the initial condi- 
tions and of the arguments to establish the results shows 
the ubiquity of thermal rectification in graded systems. 
Moreover, the existence of simple ingredients for the rec- 
tification, as described here, deserves attention: as well 
known, in the literature, the mechanism behind rectifi- 
cation in graded models is far from being clear. E.g., we 
recall the comment of G. Casati Q on the explanation 
of Chang, et al @: "the authors speculate that solitons 
might be involved in the rectification process, but this 
is still to be confirmed". Here, we do not have to make 
any speculation about the vibrational spectra or other 



2 



intricate property. 

EXISTENCE OF THERMAL RECTIFICATION 

Let us introduce our assumptions and derive our re- 
sults. 

We consider a chain with N sites, where the first site is 
connected to a thermal bath at temperature Ti , and the 
last site is connected to a bath at temperature T/v. It is 
possible to extend our analysis also for a d-dimensional 
lattice with two thermal baths at the boundaries: the 
chain structure is represented by the axis (direction) of 
the heat flow. We assume that it is possible to build 
a temperature gradient in the system. Such condition 
always happens if the Fourier's law holds, but we do not 
demand this law here (anyway, we will study cases where 
Fourier's law holds). Precisely, we assume that the heat 
flow from site j to j + 1 is given by 



j^,,+i = -/c,(vr),= 



I 



(1) 

i.e., with, say, the local thermal conductivity given by 
the average of a function of the local temperatures and 
other parameters of the system. For the homogeneous 
model, such expression reads 
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where = (l^ +Tj^+i) /2, C = 2C, which is exactly 
the formula described by several results (on homogeneous 
models) from the literature: e.g., in ref. IJ], we have 
a = 2, C'T^ = 1/K. = A^rVwV^, where A is the coefh- 
cient of the quartic anharmonic potential, uj is the coef- 
ficient of the interparticle quadratic interaction, and fj, is 
the harmonic pinning coefficient. Still for this 0^ model, 
in different conditions and methods, Bricmont and Kupi- 
ainen, [l^ and Spohn et al. [31 found ICj = T~^. And, 
in reference to works with detailed computer simulations, 
Aoki and Kusnezov [l^ obtain for this one-dimensional 
^4 model, /C cx similarly, N. Li and B. Li ob- 

tain K, (X T^^'^, with slight changes in the exponent that 
depend on the values of the pinning and anharmonic- 
ity. It is also worth to recall that, by using an analytical 
simplified scheme (derived from a rigorous and more in- 
tricate approach [19[), we obtain a similar formula for 
the local thermal conductivity of the graded anharmonic 
self-consistent chain [g^, i.e., of the anharmonic, inhomo- 
geneous model given by a chain of oscillators with quartic 
on-site potential, quadratic nearest-neighbor interparti- 
cle interaction, particles with different masses and inner 
stochastic reservoirs connected to each site. 

Let us now prove the existence of thermal rectification 
for a graded anharmonic system with a temperature gra- 
dient in the bulk, and whose local thermal conductivity 



depends on temperature (which does not follow in the 
harmonic case), and changes as we run the chain. 

From the fact that the heat current comes into the 
system by the first site, passes trough the chain and goes 
out by the last site, we have 



J'1,2 — J'2,3 — ... — J-N-l,N = J'- 
These equations together with eq.([T]) give us 

TiCiT^ + C2T.^) = T1-T2 
TiC2T^+C3T^) = T2-T3 

J-{Cn-iT^-i+CnT^) — Tn-i ~Tn. 
Summing up the equations, we find 

.(Ti-Tat) 



(2) 
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where 



K = {CiTf + 2C2T^+ . . . 

-f2CAr_iT^_i+CjvT^}~'-(iV-l), (3) 

that is the Fourier's law for the case of the thermal con- 
ductivity K, remaining finite as TV 00. From eq.([T]) and 
eq.(l2]), it follows that 



T2-T. 



CiT^ + C2T2" C2T^ + C^Tf 

— Tn 



(4) 



Thus, given the temperatures of the baths Ti and Tjv, by 
using the equations above we determine the inner tem- 
peratures T2, T3, Tn-i. For ease of computation, 
let us consider the system submitted to a small gradient 
of temperature: Ti = T + aie, TV = T + ajsie, e small. 
Hence, = T + Ofee -f 0{(!^). We wiU carry out the 
computations only up to 0{e). And so, up to 0(e), we 
have = T" + aT'^-^eak (that comes from the Taylor 
series), and 



{ak - Ofc+i) e 



CkT^ + Ck+iT^+i {Ck + Ck+i) 

as said, up to 0{e). From this equation and eq.Q, we 
obtain 



fli — 02 



a2 - a-i 



Ci +C2 C2+C3 
We may rewrite these equations as 
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Summing them up, we obtain 
ai — 02 



C1+C2 



(Ci + 2C2 



+ 2Cn-i + Cjv) — ai — ajv 



«2 = ai + j^^^ (Ci + C2) , 



where C{N) = (Ci + 2C2 + . . . + 2CAr_i + Cn)- Similarly, 
writing {ak-i~ak)/{Ck-i +Ck) instead of (ai -a2)/(Ci + 
C2) in the LHS of the list of equations above, we obtain 



ftfc = ai H = t'(fc), 

C(iV) ^ ' 



(6) 



for fc = 2, . . . , iV— 1. And so, for the thermal conductivity 
Q it follows that 

K, = [N -l)-^T°'C{N)+aT"'^e{aiCi + 2a2C2 + ■■■ 

+2aN-iCN-i+ aNCN)Y^ ■ (7) 

To investigate the existence or absence of rectification, 
we need to analyze the heat flow for the system with 
inverted thermal baths, that is, we compute the new 
thermal conductivity for the same system, but with tem- 
peratures T', where T[ = T/v and T'j^ = Ti. Following 
the previous manipulations, we see that, in the system 
with inverted baths, the new temperature for the site k 
is = T + a'l^e, where, for fc = 2, 3, . . . , - 1, 



Ofc = ajv ^7— — C(fcj. 



C{N) 



(8) 



Obviously: a[ = gn, and = oi. Hence, the expression 
for the "inverted" thermal conductivity becomes 



/C' = 



(A - 1) • |t"C(A) + aT°'-h {a[Ci + 24C2 + 



+ 2a'^_lCAr_i + a'^CAr)} . (9) 

And, with simple manipulations, we get 

1 1 _ aT"~'^e{ai ~ qn) 
IC^lC ~ (N - l)CiN) 

X {c{Nf - AQ{N) - 2CnC{N)} , (10) 

where Q{N) = C{2)C2 + ■ ■ . + C{N ~1)Cn-i- As a simple 
test for the expression above, note that it vanishes (as 
expected) in the case of a homogeneous system (Ci = 
C2 = ■ ■ ■ = Cn)- 

To continue the analysis, we take a chain with three 
sites (say, the smallest possible system). A direct com- 
putation gives us 

(C(3))2 - 4Q(3) - 2C3C(3) = Cf - Cj. 

Now we prove, by induction, that such relation is valid 
for any number of sites: we assume that it is valid for fc 



sites (i.e., for fc replacing 3 in the relation above), and 
then we show that it follows for fc -I- 1. In fact, by using 
the definitions we see that 

C(fc-t-l) = C1 + 2C2 + ... + 2Ck+Ck+i 

= C{k)+Ck+Ck+i, 
Q{k + 1) = C(2)C2 + . . . + C{k)Ck - Q(fc) + C{k)Ck. 

Then, a direct computation shows that 

{C{K + 1))2 - 4Q(fc + 1)- 2Ck+iC(k + 1) = C'l - Cl+,. 

Hence, for the difference between the thermal conductiv- 
ities of the system with N sites, we obtain 



1 1 _ aT°'-^e{ai - un) 
IC^lC ^ {N - l)C{N) 



[c!-c%], (11) 



where, we recall, aT°'^^e{ai — ajv) in the numerator 
above is T" — up to 0{e). Thus, the existence of 
thermal rectification for anisotropic, e.g. graded, mate- 
rials is transparent. 



RECTIFICATION PROPERTIES 

Now, let us examine the rectification in details and 
search for conditions leading to suitable properties. First, 
we write the expression for the rectification factor fr 



_ |/c-r| ^ iTf -T^l \Cf-C%\ 



JC' 



(C(A^))2 



Hence, fixed the temperatures at the boundaries, the be- 
havior of fr with N is given by |Cf - C%\/iCiN)f. We 
recall that 



2Cn-i + Cn « 2 / Cxdx 



N 



C{N) =Ci+ 2C2 



And, for a small gradient of temperature in the system, 

IC^iN- 1)/{T"C{N) + 0(6)}. 

Thus, to get a normal conductivity (Fourier's law) we 
must have C{N) ^ N, i.e., Cn ~ constant. That is, for 
these graded systems, at least at small temperature gra- 
dients, if the conductivity is normal then the rectification 
factor decays to zero as ^> 00. To avoid the decay of 
the rectification factor, for example, to make it finite and 
nonzero as A^ ^ 00, we need to take Cn cexp{"fN). 
And so, C(A) ~ c(exp(7A^) - For 7 > 0, C{N) has 

exponential growth and IC{N) — > as A — > 00. For 
7 < 0, C{N) constant and K. ^ N, i.e., we have 
an abnormal conductivity. That is, the regimes of non- 
decaying rectification and of normal conductivity do not 
overlap. The possibility of a non-decaying rectification is 
a very important property: as recalled before, the decay 
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of rectification is a problem for the usual diodes given by 
the sequential coupling of different parts. 

Moreover, still from the previous expression (take Ti > 
Tjv and Cn > Ci), we see that the thermal conductivity 
is smaller when the heat flows from the sites with larger 
C to the sites with smaller C. 

In short, we have shown that in a lattice system where 
it is possible to build a temperature gradient, i.e., with 
the heat flow from site j to j + 1 given by eq.(IT]), with 
graded structure (i.e. graded Cj) and with local thermal 
conductivity dependent on temperature (see eq.([T])), we 
will always have thermal rectification. To be precise, we 
need to recall that in our proof, for ease of computa- 
tion, we have assumed a system with small temperature 
gradient (however, we believe that it is not a necessary 
condition - more comments ahead). It is interesting to 
note that such conditions - temperature gradient in the 
bulk, local conductivity dependent on temperature and a 
graded structure - appear in the quantum harmonic self- 
consistent chain of oscillators 2l| , a system that presents 
rectification in opposition to its classical version (with a 
conductivity that does not depend on temperature). 

To give a concrete example, we turn to the chain 
with homogeneous anharmonic potential, homogeneous 
intcrparticle interactions, etc, but with graded masses. 
For the model with inner self-consistent reservoirs, weak 
nearest neighbor interactions, quartic anharmonicity, in 
an approximate calculation [20|, we have 



J". 



C 



where C involves the coefficients for the anharmonic- 
ity, intcrparticle interaction, etc. The denominator of 

1 /2 

the expression above may be written as [{T^ / Pj.j+i) + 

(I^Z+^/Pj+ij)], where pjj+i = TOj/toj+iTOj, pj+i,^ = 
mj+i/mj+im^-. To follow, we define 



P3 



1 (mj+i -I- mj-_i) 
2' (m,_im,+i) ' 



i.e., pj is proportional to the inverse of a reduced mass. 
Hence, considering the entire system j = 1, . . . , A^, we ap- 
proximately have J-}. j+i given by eq. (HJ withCj = 1/pjC. 
And the analysis follows as previously described: now 
with the bigger flow in the direction from the larger to 
the smaller masses. It is worth to recall that such prop- 
erty, a bigger heat flow from the larger to smaller densi- 
ties, as described here, has been already experimentally 
described 

Similar properties appear in a system with homoge- 
neous particle masses, but graded anharmonic on-site po- 
tentials or graded interpaticle interactions. 



FINAL REMARKS 

We have some remarks. First, we stress that we have 
presented here sufficient, not necessary, conditions for 
manifesting thermal rectification in anisotropic systems. 
In ref.[22j|. by computer simulations, the authors describe 
rectification in a graded mass Fermi-Pasta-Ulam chain, a 
model with an invariant translational potential and ab- 
normal conductivity (even for the case of homogenous 
masses). We also recall that, for the (very different) case 
of a system of two-terminal junctions, sufficient condi- 
tions for rectification have been described in a recent 



work by Wu and Segal [23 1. 

A further investigation of great interest is the behav- 
ior of the graded system as submitted to a large gradient 
of temperature: we believe that it shall lead to a signi- 
ficative rectification. In ref. 2J] , for some specific models 
given by chaotic billiard systems, the authors claim that 
there is a significative rectification "provided the tem- 
peratures (of the two sides of the system) are strongly 
different..." 

To conclude, we emphasize that due to their simplicity, 
the assumptions and arguments described here follow for 
many of the usual systems modeling heat conduction in 
solids: it shows the ubiquity of rectification in graded 
systems. Moreover, the existence of simple conditions for 
the existence of an efficient rectification, and the fact that 
graded systems may be constructed in practice (and are 
even abundant in nature) indicate that they are optimal 
material to be used in the construction of a thermal diode 
(and also thermal transistors, etc), and so, their use shall 
certainly contribute to the advance of phononics. 

Work supported by CNPq (Brazil). 
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